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1.  Introduction 


In  a recent  investigation  a system  of  linear  elastic  equations  for  a 
two-constituent  composite  material  was  obtained  from  a general  nonlinear  system 
for  N-constituents.  The  general  system  of  nonlinear  equations  was  obtained 
from  a model  of  the  composite  consisting  of  interpenetrating  solid  continua, 
in  which  the  motion  of  a point  of  the  combined  continuum  could  be  finite  while 
the  relative  motion  of  each  of  the  constituents  was  constrained  to  be  in- 
finitesimal in  order  that  the  solid  composite  not  rupture.  The  aforementioned 
linear  equations  were  written  in  explicit  detail  for  the  isotropic  and  trans- 
versely isotropic  symmetries.  The  linear  elastic  equations  for  the  two- 
constituent  transversely  isotropic  composite  obtained  in  the  earlier  work1 
form  the  basis  of  the  work  presented  here. 

In  this  paper  the  aforementioned  linear  equations  for  the  transversely 
isotropic  composite  material  are  applied  in  the  analysis  of  some  simple  but 
very  interesting  one-dimensional  static  load  transfer  problems  and  the  propa- 
gation of  straight-crested  surface  waves.  In  the  two  particular  one- 
dimensional static  problems  considered, axial  loading  is  applied  to  the  fiber 
reinforcement  which  enters  the  unloaded  matrix.  In  one  case  both  the  matrix 
and  reinforcement  are  held  fixed  at  the  supporting  end,  while  in  the  other 
case  the  reinforcement  ends  at  some  distance  into  the  matrix.  In  both 
cases  the  influence  of  gravity  is  included  in  the  analysis.  In  each 
instance  the  rate  of  transfer  of  stress  from  the  reinforcement  to  the  matrix 
is  determined  in  terms  of  the  defined  material  coefficients  of  the  two- 
constituent  composite  as  a byproduct  of  the  solution  of  a simple  system  of 
ordinary  differential  equations  with  constant  coefficients.  The  treatment 
of  such  problems  within  the  framework  of  the  theory  of  linear  elasticity  is 

prohibitively  complicated  and  such  problems  cannot  even  be  mathematically 
defined  using  the  ordinary  strength  of  materials  approach. 


In  the  case  of  surface  wave  propagation  only  the  placement  of  fiber 

reinforcement  parallel  to  the  free  surface  is  considered.  The  formal  solution 

for  surface  wave  propagation  in  the  direction  of  the  fiber  reinforcement  is 

presented.  However,  since  the  material  coefficients  occurring  in  the  theory 

have  never  been  measured  for  any  composite,  a calculation  cannot  be  performed. 

Nevertheless,  if  the  resulting  theory  is  reduced  by  making  the  plausible 

simplifying  physical  assumptions,  for  certain  cases  of  interest,  that  only 

axial  stress  exists  in  the  reinforcement  and  the  entire  interaction  depends 

only  on  the  relative  displacements  of  the  constituents  in  the  model,  the 

remaining  unknown  constants  in  the  description  can  be  partially  estimated  from 

the  known  ordinary  elastic  constants.  Then  by  making  an  additional  assumption, 

calculations  can  be  performed.  At  this  point  it  should  be  noted  that  when  the 

aforementioned  simplification  is  made,  the  resulting  linear  equations  are 

2 3 

identical  with  an  earlier  system  due  to  Bedford  and  Stern  ’ . The  aforemen- 
tioned partial  estimation  of  the  unknown  constants  in  the  simplified  descrip- 

4 

tion  is  made  using  a procedure  due  to  Martin,  Bedford  and  Stern  . With  the 
material  constants  thus  determined  the  dispersive  surface  wave  velocity  has 
been  calculated  for  a glass  fiber  reinforced  phenolic  resin.  In  this  simpli- 
fied theory  upper  (in  frequency)  optical  type  surface  wave  branches  are  found 
in  addition  to  the  lower  acoustic  type  surface  wave  branches.  It  should  be 
mentioned  that  it  is  felt  that  the  upper  surface  wave  branch  that  occurs  in 
the  simplified  description  will  not  occur  in  the  full  description  because  in 
the  latter  case  all  the  independent  solutions  of  the  differential  equations 
that  remain  coupled  by  the  boundary  conditions  will  probably  not  decay  with 
depth.  In  each  case  treated  the  acoustic  type  surface  wave  branch  turns  out 
to  be  asymptotic  to  the  non-dispersive  surface  wave  velocity  of  the  matrix  at 


very  long  wavelengths.  It  should  be  noted  that  the  theory  employed  and,  of 


3. 


course,  solutions  presented  are  valid  only  for  wavelengths  long  compared  to 
the  spacing  of  the  fiber  reinforcement  in  much  the  same  manner  that  the  theory 
of  elasticity  is  valid  only  for  wavelengths  large  compared  to  a lattice 
spacing.  However,  because  of  the  behavior  of  the  branches,  we  obtain  and 
present  results  considerably  beyond  the  range  of  validity  of  the  theory. 
Nevertheless,  we  indicate  the  limit  of  validity  of  the  theory  on  each  branch 
plotted. 


2.  One-Dimensional  Static  Problems 

In  this  section  we  consider  some  simple  but  interesting  one-dimensional 
static  problems  for  two-constituent  transversely  isotropic  composite  materials. 
In  each  instance  the  load  is  applied  in  the  preferred  direction  of  transverse 
isotropy,  which  lies  along  the  length  of  the  parallel  fibers,  and  it  is 
assumed  that  all  displacement  and  relative  displacement  components  transverse 
to  this  direction  are  constrained  to  vanish  and  that  the  remaining  displace- 
ment variables  are  independent  of  the  transverse  coordinates.  We  consider 
only  problems  that  satisfy  these  criteria.  Under  these  circumstances  the 
nontrivial  linear  constitutive  equations  take  the  form3 

* * ( 1 ) 

>il  * K22  * V3,3  + *3*3,3-  ,2'1> 

. A (1  \ 

K33  = C5U3,3  + P6W3,3’  (2,2) 

„ A - (1) 

^11  = ^22  = P5U3,3 +b3W3,3  ’ (2,3) 

■*33  = ^6U3,3  + b5W3,  3 ’ (2,4) 

3 = - a2w3(1)  , (2.5) 

where  is  the  preferred  direction  of  transverse  isotropy,  u3  is  the  non- 
zero displacement  component  of  the  center  of  mass  of  the  combined  two- 
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constituent  composite  material  and  w^^  is  non-zero  component  of  the 

relative  displacement  of  the  continuum  representing  the  matrix.  We  employ 

Cartesian  tensor  notation  and  we  have  introduced  the  convention  that  a comma 

followed  by  an  index  denotes  partial  differentiation  with  respect  to  a space 

coordinate.  The  represent  the  components  of  the  stress  tensor  for  the 

combined  continuum  and  the  3 „ represent  the  relative  stress  tensor  which  is 

LM 


defined  by 


P=P 


<1,  +pC>( 


(2.6) 


where  T ^ and  p ^ represent  the  components  of  the  stress  tensor  and  mass 
LM 

density,  respectively,  of  each  of  the  interpenetrating  continua.  The  vector 

field  3L  is  related  to  the  volumetric  force  of  interaction  between  the  two 
M 


constituents  by  the  relation 


3M  ■ LpM2  <l  + '»  ’ 


(2.7) 


L 12 

where  F„  is  the  volumetric  force  exerted  by  continuum  2 on  continuum  1.  At 
M 

this  point  it  is  to  be  noted  that  the  T and  p ^ do  not  represent  the  actual 
components  of  stress  and  mass  density  of  each  of  the  constituents  in  the 
composite,  but  only  represent  those  quantities  in  each  of  the  interpenetrating 
continua,  which  occupy  the  same  region  of  space  and,  respectively,  represent 
each  constituent  in  the  model.  As  a consequence,  if  A™  and  A^  represent  the 
areas  occupied  by  the  matrix  and  fibers,  respectively,  in  a typical  area  A 
of  the  interpenetrating  continua  normal  to  the  fiber  length,  we  have 


m , f (1)  mm.  (2)  f f „ 

A = A + A , p ' ' = p A /A  , p = p A /A  , 


m 
r . 


. =T.(1)A/Am,  Tf.  =T.(2)A/Af, 


il 


il 


(2.8) 


where  the  variables  with  the  superscripts  m and  f represent  the  actual  re- 
spective quantities  in  the  matrix  and  fiber  reinforcement,  respectively. 


5. 


The  remaining  nontrivial  stress  equations  of  equilibrium  and  relative  stress 


equations  of  equilibrium  are 


*33.3  + Pf3  = 


j)  + 5 + p ' ;f  = 0 


(2.10) 


where 


-pV+pV,  -f,(2)  , 


(2.11) 


and  and  f.^  denote  the  components  of  body  force  per  unit  mass  in  the 

continua  representing  the  matrix  and  fiber  reinforcement,  respectively,  which 
in  the  case  of  the  gravity  force  are  the  same  as  the  body  force  intensities  f” 
in  the  matrix  and  f^  in  the  fiber  reinforcement,  both  of  which  equal  g. 


The  substitution  of  (2.2),  (2.4),  (2.5)  and  (2.11)  into  (2.9)  and  (2.10) 


yields 


C5U3,33  + P6W3,33  + pg  = ° ’ 
P6U3,33  + b5W3.33  “ a2W3  = °’ 


(2.12) 


(2.13) 


which  are  the  one-dimensional  displacement  equations  of  equilibrium  that  apply 
to  the  one-dimensional  static  problems  treated  in  this  section.  The  solution 
to  (2.12)  and  (2.13)  may  be  written  in  the  form 


P6  -®x3  0*3  i 2 

u = - — (Ae  + De  ) + pgx  + Bx  + C , 

c5  2c  5 

(1)  _Q(X3  “*3 

w ' ' = Ae  + De  + yg  , 


(2.14) 


2 /\  A . A A i 4 n A 

Q = c5a2/(c5b5  - P6)  , V=Pe6/c5a2, 


(2.15) 


and  A,  B,  C and  D are  arbitrary  constants,  to  be  found  by  satisfying  boundary 
conditions  in  a given  one-dimensional  problem. 


6. 


In  the  first  specific  problem  we  consider,  a total  compressive  force  P is 
applied  to  all  the  fibers  crossing  a given  cross-sectional  area  of  a combined 
composite  material  body  of  height  i resting  on  a rigid  surface,  as  shown  in 
Fig.l.  Since  the  supporting  surface  is  rigid,  there  is  no  displacement  in 
either  the  matrix  or  fibers  at  the  support  and  we  have 


u3  + w3(1)  =0,  u3  + wf>  = 0 , at  x3  = 0 . (2.16) 

Since  u^  is  defined  as  the  displacement  of  the  center  of  mass  of  a point  of 
the  combined  composite  continuum,  we  have 


CD  (D 

0 w3 


_ (2)  (2) 
P 


= 0 


(2.17) 


which  with  (2.16)  enables  us  to  write  the  boundary  conditions 


U3  = 0 , w<X>  = 0 , at  x3  = 0 . 


(2.18) 


Since  no  force  is  applied  to  the  matrix  at  x3  = £,  we  have 


Td)  _ o T<2)  =__  at  x , l 

T33  " ° ’ 33  Po  ’ 3 Z> 


(2.19) 


where  p = p/a  and  which  with  (2.6)  and 
o 

(1)  (2) 

K = T + T 
L.M  LM  LM 


(2.20) 


enables  us  to  write 


K33  = ~ Po  ’ ^33  = rpo  ' at  x3  = £ * 


Now,  the  substitution  of  (2.14)  into  (2.18)  and  (2.21)  yields 


A = - 


Yg  + Kp^e 


-ai 


-2  oil  ' 


D = 


Yg  - up  e 
2a  l 


aH 


1 + e 

B = - (pg^+Pc)/c5,  C = - ($6/c5)\g, 


1 + e 


where 


(2.21) 


(2. 22) 


h = «*/a2) (r  + P6/c5)  . 


(2.23) 


7. 


The  substitution  of  (2.22),  with  (2.23),  in  (2.14)  yields  the  solution.  The 

substitution  of  (2. 14), with  (2.22)  and  (2.23),  into  (2.1)  - (2.5)  yields  all 

the  stresses,  relative  stresses  and  interaction  forces,  which  we  do  not  bother 

to  write.  However,  in  the  absence  of  g for  the  actual  stresses  in  the 

fibers  and  t™  in  the  matrix  at  the  support  = 0,  we  obtain 

2 

m _ ^0  A T ^6.  6.  r \ 1 

T33  1 + r m Lr  + cr  + \c_  5 cosh  aHj  ’ 


1 

0 

a r ^6 

(h 

1 I V- 

— r + — 

+ a\T 

-L  t r 

Am  L C5 

VC5 

1 

0 

a r.  h. 

%2 

-Qfe 

1 + r 

Af  L C5 

Vc5 

) — 

>/  cos! 


5/  cosh  all 


(2.24) 


In  the  second  specific  problem. we  consider  fiber  reinforcement  entering 
a matrix  and  terminating  uniformly  at  a distance  i into  the  matrix  which  con- 
tinues down  to  a rigid  support  at  a distance  b below  the  junction,  as  shown 
in  Fig. 2.  A total  tensile  force  P is  applied  to  all  the  fibers  crossing  a 
given  cross-sectional  area.  Since  the  continuum  representing  the  matrix  and 
the  continuum  representing  the  fibers  can  neither  separate  from  nor  penetrate 
into  the  single  matrix  continuum  that  abuts  the  composite  at  the  junction, 
we  must  have 


w3  = 0 , w3  = 0 , at  x3  = 0 , 


(2.25) 


which  is  consistent  with  (2.17).  In  addition,  the  displacement  u3  of  the 
center  of  mass  of  the  combined  composite  continuum  must  be  the  same  as  the 
displacement  U3  of  the  isotropic  single  matrix  continuum  at  the  junction. 
Consequently,  as  kinematic  boundary  conditions  at  the  junction  we  have 


where  U3  satisfies 


w3  = 0 , u3  = U3  , at  x3  = 0 , 


<\m  + 2nm)u3  33  - pmg  = 0 . 


(2.26) 


(2. 27) 


8. 


and  the  non-trivial  stress  components  in  the  single  matrix  continuum  are  given 
by 


mm  ,,m  . _ m.  m m , m 

33  - ( 2p.  )u3;3  * tl1  ~ t22  _ X U3,3  * 


(2. 28) 


In  addition  to  the  continuity  of  displacement  at  = 0 we  have  the  continuity 


of  traction,  i.e., 

*33  = T33  ’ at  X3  = 0 ' 

Since  no  force  is  applied  to  the  matrix  at  x3  = lt  we  have 

„ (1) 


■33  = ° » T33  = Po  ’ at  x3  = 


(2)  _ 

33 

which  with  (2.6)  enables  us  to  write  the  boundary  conditions 


K33  = Po>  ^33  = - rPo>  at  X3  = l’ 
Since  the  supporting  surface  is  rigid,  we  have 


U3  = 0 , at  x3  = - b . 


(2. 29) 


(2.30) 


(2.31) 


(2.32) 


Thus,  the  boundary  conditions  are  (2.26),  (2.29),  (2.31)  and  (2.32).  The 

solution  to  (2.27)  may  be  written  in  the  form 


U3  = [pmg/2(Xm  + 2Mm)]  x3  + Ex3  + F. 


(2.33) 


Now,  the  substitution  of  (2.14)  and  (2.33)  into  (2.26),  (2.29),  (2.31)  and 


(2.32)  yields 


A = - 


Yg  + np  e 
o 


-a  l 


1 + e 


-2  0/1 


D = - 


Yg  - xPQe 


a l 


1 + e 


2 a l 


i (P  - pg-2  - 4 Pmgb)b 

B = 4 (p  - pg£)  , c = - -|  Yg  + —2 ± 

c c X + 2m 


E = 


(po  - pgi) 
X + 2m 


F = 


(pQ  - pg^  - Pmgb)b 
. m _ m 

X + 2m 


(2.34) 


which  when  substituted  in  (2.14)  and  (2.33),  respectively,  yields  the  solu- 
tion. The  substitution  of  (2.14)  and  (2.33)  with  (2.34)  into  (2.1)-  (2.5) 


.....tt,!  ..LJ.,  HIP  WIW  ^ 


9. 


and  (2.28) , respectively,  yields  all  the  stresses,  relative  stresses  and 
interaction  forces,  which  we  do  not  bother  to  write.  However,  in  the  absence 
of  g for  the  actual  stresses  in  the  fibers  and  in  the  matrix  at  the 
junction  at  x^  = 0,  we  obtain 


m 

Po 

A 

t-  i - 4s 

- — ) 

H 

—1 

33 

1 + r 

m 

A 

cs 

cosh 

f 

Po 

A 

t1  - k + 4s 
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H 

-i 

33 

1 + r 

Af 

V 

cosh 

oil _ 

(2.35) 


Using  a highly  simplified  model,  the  results  of  this  analysis  are  presented 
in  the  Appendix  in  terms  of  the  known  material  constants  of  the  fiber  rein- 
forcement and  matrix  and  the  relative  geometry. 


3 . Surface  Waves 

In  this  section  we  consider  surface  waves  propagating  along  the  free- 
surface  of  a semi-inf inite  fiber  reinforced  composite  material.  Only  the 
case  of  straight-crested  surface  waves  propagating  in  the  direction  of  the 
fiber  reinforcement,  which  runs  parallel  to  the  free-surface,  is  considered 
because  the  method  of  analysis  is  essentially  the  same  for  surface  waves 
propagating  normal  to  the  direction  of  the  fiber  reinforcement,  or  in  any 
other  direction  for  that  matter,  and  calculations  based  on  the  analysis  are 
not  performed  because  the  material  constants  are  not  presently  known  for  any 
two-constituent  composite  material.  However,  in  the  next  section  calculations 
are  performed  for  surface  waves  propagating  both  in  and  normal  to  the  direction 
of  the  fiber  reinforcement  using  a highly  simplified  model  of  the  material, 
and  it  is  indicated  that  certain  of  the  surface  waves  allowed  by  the  simplified 
model  will  not  be  permitted  by  the  more  general  model  considered  in  this  section. 


10. 


A schematic  diagram  of  the  free-surface  of  the  two-constituent  composite 
along  with  the  associated  coordinate  system  is  shown  in  Fig. 3.  Since  we 
consider  straight-crested  surface  waves  propagating  in  the  direction  of  the 
fiber  reinforcement,  which  is  in  the  x^-direction,  the  solution  functions 
are  independent  of  x^.  Furthermore,  an  examination  of  the  equations  for  the 
two-constituent  transversely  isotropic  composite  with  x^  the  preferred 
direction5  presented  in  Ref.l  reveals  that  for  x^~ independence  u 2 and  w 
uncouple  from  u^,  w , u^  and  w^.  Hence,  under  these  circumstances  u?  and 
may  be  taken  to  vanish  and  the  non-trivial  differential  equations  may  be 
written  in  the  form 

(C1  + 2C3)U1,11  + C4U1,33  + (P1  + 32)Wi;LU  + 1 P4W1^33  + 


(C2  + C4)U3,13  + (I  P4  + = P”l  ’ f3‘ 

31  + e2)Ul,ll  + 1 P4U1,33  + (2bl  + b2)Wl!ll  + Vb4  + 1 b7X(133  + 

(I  e4  + 0U3,13  + (b3  + b4  " I b7>M3“alW?)  -rpW1>  ’ (3 ' 

(C2  + C4)U1,13  + (l  P4  + + C4U3,11  + C5U3,33  + 


(c2  + C4}' 


13  + C4U3,11  + C5U3,33  + 


2 P4W3,11  + P6W3,33  pU3  > 

(i  04  + P3K,13  + (b3  + b4  - \ b7>M3  + \ V3,ll  + 
P6U3,33  + (b4  + 2 b7)W3,ll+b5W3(,33_a2W3  = rpW3  * » 

and  the  non-trivial  boundary  conditions  may  be  written  in  the  form 

hi  * hi  "hi  ‘ *13  - 0 “ *1  ’ 

where 

*11=  (c1  + 2c3)u11+  (P1+P2)w^+  (c1  + c2)u3;3+  (g1  + P3)w3a>, 


^11  " (P1  + P2)U1,1+  (2bl  +b2)Wl,  !+  (P1  + P5)U3,3+  (b2+b3)W3,3  > 
K13  " C4(U3,1+U1,3)+  1 P4C,,3(^+WM)» 

^13  = 2 g4  (ul,3  +U3,1)  + (b4  " 2 ^^1,1  + lb4  + 


and  everything  vanishes  as  x^  -*  °°. 

As  a solution  of  the  differential  equations  we  take 


i(T]x1+§x3-u)t) 

u<*  = V ’ wa  = 32e 


i (T|x  +§x3-u)t) 


, or=l,3,  (3.10) 


which  satisfies  (3.1)  - (3.4)  provided 


®11  - >\  + C12E1  + =13*3  * 514B3  * ° > 


a21Al  + <«22  - lBl  * ®23S3  + 824®3  ' ° • 


«3lh  + 232B1  + ®33  - >”  ,fl3  + °34B3-°  • 
®4lh  * 542B1  * °43a3  + ®44  ‘ 


where 


Qll  = (Cl  + 2c3)T12  + c4?2 , q12  = q21  = (P,  t P2)T12  + i P452  , 

213  = 231  = (C2  + C4)T1^  ’ 214  = 241  (l  ^4  + ’ 

222  (2b1+b2)Tl  + (b4  + 7 b7)5  + ai  > 223  = C32  = (l  P4  + 0^  ’ 

Q24  = Q42  = (b3+b4  - khl)^>  C33  = ^ + V1"  ’ 

^34  - *43  = I + > 244  - (b4  + ± b?K  ^ V2  - a2  . 


Equations  (3.11)  constitute  a system  of  four  linear  homogeneous  algebraic 
equations  in  B^,  A3  and  B3,  which  yields  non- trivial  solutions  when  the 
determinant  of  the  coefficients  of  A^,  B^,  A3  and  B3  vanishes,  i.e.,  when 
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2 

- pU>  ) 

2l2  2 

213 

214 

21 

(Q22  ~ rPUJ  ) 

223 

C24 

= 0 . 

(3 

31 

232 

(Q33  - pco  ) 

Q34  2 
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Q42 
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(Q44  - rpu>  ) 

2 2 2 

Equation  (3.13)  is  a quartic  in  a)  , ? and  T|  and,  hence,  for  a given  cu  and  § 

2 

there  are  four  in  general  complex  T]  . Since  the  solution  functions  must  vanish 
as  x^  -»  only  those  T|  ^ (n  = 1,  2,3,4)  with  positive  imaginary  part  are 
admissible.  For  a given  T]  , three  of  the  four  equations  in  (3.11)  yield 
amplitude  ratios,  which  we  denote 


A™  : B]L(n)  : A3(n)  , sf* 


(3.14) 


All  four  admissible  solutions  at  a given  ^ and  UJ  are  required  in  order  to 
satisfy  the  four  boundary  conditions  in  (3.5).  Consequently,  we  take 

r (n)  (n)  iT1<n)xl 
u = / C A e e 


a 


n=l 


a 


(1)  _ f c (n)  (n)  i1'l”>xl  ‘'V11 

wq  " L C Ba  e e 

n=l 


(3.15) 


which  satisfy  (3.5)  provided 
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; C(n)L^n)  = 0 , Y = 1,  2,3,4  , 
n=l 


(3.16) 


where 
(n) 


!•£"'  = (c1  + 2c3)Tl(n)A1(n)  + (P1  + P2)Tl(n)B1(n)  + (C1  + c2)?A3(n)  + (^  + ^ 

L2n)  = (P1  + P2)'n(n)A^n>  + (2b3  +b2)T](n)B1(n)  + (@1+ 35)5A3(n)  + (b2  + b3)5B3(n) 

t Cn)  - r Fa  + 1 R Fn  ( _ -.(n),  (n)  1 R P (n)  (n) 

L3  - C4§A1  + 2 345B1  + c4*  A3  + 2 W B3  ’ 

T (*"*)_  ^*  q crn  ^ ^ \ eD  (^ ) , 1 o ri  (t  )fl(n)  1 . .„(n)  (n) 

L4  “2  345A1  +(b4_2b7)5Bl  + 2 M A3  + (b4  + 2b7)T>  B3  * 


(3.17) 


k 


) 

Equations  (3.16)  constitute  four  linear  homogeneous  algebraic  equations  in 
the  and  for  a non-trivial  solution  the  determinant  of  the  coefficients 

must  vanish,  i.e., 


(1) 

1 

h<3’ 

(1) 

'2 

(1) 

'3 

(1) 
' A 

l.(2> 

i.;31 

i.!4’ 

If  the  constants  are  known,  a calculation  proceeds  by  selecting  values 
for  § and  uu,  which  enables  the  determination  of  the  four  T)^  from  (3.13)  and 
the  attendant  amplitude  ratios  from  (3.11),  Then  everything  in  (3.18)  is  known 
and  either  it  is  satisfied  or  it  is  not.  If  (3.18)  is  not  satisfied,  change 
either  § or  uu  and  repeat  the  entire  calculation  until  (3.18)  is  satisfied. 

Experience  with  surface  waves  indicates  that  there  will  be  one  cu  at  a given  c. 
that  satisfies  all  conditions,  i.e.,  (3.13)  and  (3.18),  However,  since  the 

constants  are  not  known  for  any  two-constituent  composite  material,  a calcula- 
tion cannot  be  performed.  In  the  next  section  calculations  are  performed  for 
a highly  simplified  version  of  the  model  for  which  the  constants  can  be  esti- 
mated from  the  known  constants  of  the  two  constituents  in  the  composite. 

4 . Surface  Waves  and  the  Simplified  Model 

In  this  section  we  simplify  the  equations  for  the  two-constituent  trans- 
versely isotropic  composite  by  reducing  the  model  sufficiently  that  the 
material  constants  of  the  composite  can  be  estimated  from  the  known  constants 
of  the  individual  constituents  of  the  composite  while  still  retaining  certain 
of  the  essential  characteristics  of  the  composite.  Since  in  the  simplified 
model  the  constants  are  essentially  known,  calculations  can  be  and,  indeed, 
are  performed  for  surface  waves  propagating  in  two-constituent  transversely 

1 


14. 


isotropic  composite  materials.  The  simplified  model  we  are  considering  is 
for  a fiber  reinforced  composite  material  consisting  of  an  elastic  matrix 
containing  uniformly  distributed  continuous  fibers  extending  in  the  x^-direction, 
in  which  the  fibers  occupy  a small  fraction  of  the  total  composite  volume. 

On  account  of  the  latter  condition  in  the  reduced  model  it  is  assumed  that  the 
stresses  in  the  matrix  are  related  to  the  strains  in  the  matrix  by  the  consti- 
tutive relations  of  linear  isotropic  elasticity  and  are  independent  of  the 
strains  in  the  fibers.  Similarly,  the  stresses  in  the  fibers  are  assumed  to 

be  independent  of  the  strains  in  the  matrix.  It  is  further  assumed  that  all 

(2) 

stress  components  in  the  long  narrow  fibers  vanish  save  the  axial  stress  . 

which  may  then  be  written  as  a function  of  the  axial  strain  in  the  fibers  only. 

Although  this  latter  assumption  seems  questionable  to  us  in  general  we  make  it 

anyway.  Then  the  only  interaction  between  the  matrix  and  the  fibers  remaining 

L 12 

is  the  volumetric  interaction  term  F which  from  (2.5)  and  (2.7)  takes  the 

M 


LFp2  =-  (1  + t)'1a1wp1)  , L?32  =-  (1  + r)_1a2w3(1)  , P = 1,2.  (4.1) 

At  this  point  it  should  be  noted  that  the  aforementioned  assumptions  make  the 

reduced  model  for  the  linear  case  identical  with  that  of  Martin,  Bedford  and 
4 

Stern  . 

On  account  of  the  assumptions  made  in  the  simplified  model  of  the  two- 

constituent  composite,  it  is  advantageous  to  write  the  equations  in  terms  of 

the  infinitesimal  displacement  fields  of  each  constituent  instead  of  the 

£ 

center  of  mass  and  relative  displacement  fields.  To  this  end  we  write 


(1)  (1)  (2)  (2) 

U K = UK  + WK  > UK  = UK  + WK  • 


which  with  the  three  -dimensional  version  of  (2.17) 


15. 


(D  - w(2) 

™K  -"WK  ’ 


enables  us  to  write 


-1 


uK  = (1  + r)  (ruR 


(1) 


(2) 


.(1) 


-1,..  (1) 


+ uK  ) , wK  = (1  + r)  (u 


- u, 


(2) 


(4.3) 


(4.4) 


From  (2.6)^  and  (2.20)  we  obtain 


(1) 


(4.5) 


Substituting  from  (4.4)  into  the  constitutive  equations  for  the  two-constituent 
7 

isotropic  composite  and  then  into  (4.5).  and  introducing  the  aforementioned 
assumptions  of  independence  of  the  stress  in  the  matrix  continuum  on  the  strain 
in  the  fiber  continuum  and  vice-versa,  which  introduces  relations  among  the 
more  generally  defined  material  constants,  and  the  further  assumption  of  uni- 
axial stress  in  the  fiber  continuum,  we  obtain 


t(2>  -e,2,u(2)5  . 

LM  3,3  LM  ’ 


(4.6) 


where  X ^ and  p ^ ^ are  the  Lame  constants  of  the  matrix  continuum  and  is 

Young's  modulus  for  the  fiber  continuum,  which  are  related  to  the  respective 
constants  in  the  matrix  and  fibers  by 

X{1)  = AV/A  , p(1)  =pmAm/A,  E(2)  =EfAf/A.  (4.7) 

Since  the  fiber  reinforcement  in  this  simplified  model  is  restricted  to  occupy 

a small  fraction  of  the  total  composite  volume,  (4.7)  enables  us  to  write 

i (D  (1)  m „ (2)  „f  f ,,  , 

X = X , u.  = |1  , E = E NS  , (4.8) 

where  S^  is  the  cross-sectional  area  of  each  fiber  and  N is  the  number  of 
fibers  per  unit  area.  In  order  to  complete  the  constitutive  equations  for  the 


reduced  model  we  substitute  from  (4.4)  into  (4.1)  to  obtain 


L 12  _ ~ , (1)  (2) 

F3  “ a2  U3  " U3  5 


where 


= a^/  (1  + r)2,  a2  = a2/ (1  + r)2  . 


(4.10) 


"The  stress  and  relative  stress  equations  of  motion  in  the  absence  of 
body  and  relative  body  forces  take  the  respective  forms 


Vl  = PV  ^LM,L  + = rf% 


(4.11) 


which  with  (2.6),  (2.7),  (2.20)  and  (4.4)  enables  us  to  write 


T<1>  + V2  = pa)u{1) 

LM,  L M K M 

T<2>  - V2  = p(2)u(2) 

LM,  L M M 


(4.12) 


(4.13) 


the  latter  of  which,  on  account  of  the  assumption  of  uniaxial  stress  in  the 


fibers  yields 


(2)  L 12  (2)  ••  (2) 

T33,36LM"  FM  = P UM 


(4.14) 


Substituting  from  (4.6)  and  (4.9)  into  (4.12)  and  (4.14)  and  employing  (4.8), 


we  obtain 


m m,  (1)  m ~ , (2)  (1), 

(X  +u)u'  +uu  + a (u  -u'  ) 

V ^ ' K.KL  P L,KK  1V  L L 

~ w,(2)  (1 ) . *>  _ „(1)..(D 

+ (a2  - al)(uL  - uL  )63L-p  ul  . 


(4.15) 


which  are  the  displacement  equations  of  motion  of  this  highly  simplified  model 

of  the  two-constituent  composite  material.  Equations  (4.15)  and  (4.16)  are 

4 

identical  with  the  equations  of  Martin,  Bedford  and  Stern  , who  have  provided 


an  approximate  procedure  for  estimating  the  constant  a2  in  terms  of  the  known 
constants  of  the  individual  constituents  in  the  composite  and  the  geometry. 


17. 


Their  analysis  provides  the  result 


a2  = 


(/Vh2)  (1  - 62/h2)2 


(4. 17) 


i log  (h/6 ) +~  62/h2 


1 r4  .4  3~|  1 m f .2  2 2 

--6/h  - - +-(p,/p,)(l-6/h) 


where 

h2  = (/3/2tt ) s2  , 


(4. 18) 


and  6 is  the  fiber  radius,  h is  the  radius  of  a cylinder,  each  of  which  encloses 
a single  fiber  in  the  hexagonal  array  and  abuts  all  adjacent  cylinders,  s is  the 
fiber  spacing  and  p,™  and  |i^  are  the  shear  moduli  of  the  matrix  and  fiber 
material,  respectively.  However,  the  constant  a^  still  remains  undetermined 
in  this  simplified  description  and  it  does  not  appear  to  be  possible  to 
estimate  it  in  some  approximate  manner  because  an  appropriate  problem  yielding 
a simple  solution  of  a full  linearly  elastic  boundary  value  problem  cannot  be 
found.  Nevertheless,  since  we  are  introducing  this  highly  limited  simplified 
description  in  order  to  obtain  some  numerical  results,  we  take  the  rather 
arbitrary  course  of  assuming  that  a^  = a2.  Under  these  circumstances  all 
material  constants  in  the  equations  are  known  and  calculations  can  readily  be 
performed. 

As  in  the  case  of  the  more  general  model,  in  the  consideration  of  straight- 
crested  surface  waves  propagating  in  the  direction  of  the  fiber  reinforcement 
^-direction),  the  solution  functions  are  independent  of  x^  and  u^1^  and  uj2^ 
vanish.  Under  these  circumstances  we  have  the  differential  equations  (4.15) 
and  (4.16),  with  x - independence  understood  and  u^  * = u.|2^  = ai°n9  with 
the  non-trivial  boundary  conditions 


T 


(1) 

11 


at  = 0 . 


(4. 19) 


(2)  (2) 

Since  and  t vanish  identically  in  the  simplified  model,  we  have  only 

the  two  non-trivial  boundary  conditions 


18. 


- 'S’  - » 


at  Xj  = 0 . 


(4.20) 


It  is  our  belief  that  the  foregoing  is  essentially  the  reason  that  the  surface 
wave  solutions  in  the  simplified  model  yield  upper  (optical  type)  surface  wave 
branches,  which  we  feel  will  not  be  in  conformity  with  results  that  would  be 
obtained  from  the  more  general  model. 

As  a solution  of  the  differential  equations  we  take 


(i)  _ 1 (Tlx1+?x3-U)t)  (2)  i OlXj+^-uut) 


u = A e 
O'  O' 


u =»  B e 
> a a e 


, *=1,3, 


which  satisfies  (4.15)  and  (4.16)  provided 


(4.21) 


where 


(1)  2 ~ 

(P1L  - p U)  )A1  - a1B1  + P13A3  = 0, 

a/  / 2 ) 2 

- \)A1  + ^ - p 'cu  )B1  = 0 , 

(1)  2 ~ 

P31A1  + (P33  - P « )A3  - a2B3  = °> 

~ (2)2 
‘ V3  + (P44  " P <"  )B3  = °- 


pn  = (Xm+2nm)Tl2  + nm|2  + al,  P13  = P31=  {\m  + n")T|5  , 


m_2 


,m  „ m, _2 


P = + 4 T]  + (X  +24  )r  + a„  P,,=Ev*'r  + 


= .(2)P2  A ~ 


(4.22) 


(4.23) 


33  r-  , “2'  ‘44  “ J “2 

Equations  (4.22)  constitute  a system  of  linear,  homogeneous  algebraic  equations 
in  Ap  B^,  A3  and  B3,  which  yields  non-trivial  solutions  when  the  determinant 
of  the  coefficients  of  Ap  Ap  B^  and  B3  vanishes,  i.e.,  when 


2 2 2 

Equation  (4.24)  is  quadratic  in  T)  , cubic  in  § and  quartic  in  cu  , and,  hence, 

2 

for  a given  uu  and  § there  are  two  in  general  complex  T]  . Since  the  solution 

/ „ \ 

functions  must  vanish  as  x^^  -*  only  those  T]  ' (n=l,  2)  with  positive 

/_  \ 

imaginary  part  are  admissible.  For  a given  T1  three  of  the  four  equations 
in  (4.22)  yield  amplitude  ratios,  which  we  denote 


A(n)  • B(n)  • A(n)  • B(n) 
1 1 • A3  • 3 


Since  only  two  non-trivial  boundary  conditions  remain  in  (4.20),  the  two 


admissible  solutions  are  adequate  and  we  take 

(1)  V i(5x3-“t> 

u = / C A e e 

o (La 

n=l 
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At  this  point  it  should  be  noted  that  on  account  of  the  simplified  model 
Eqs.  (4.27)  do  not  contain  the  . Equations  (4.27)  constitute  two  linear 

homogeneous  algebraic  equations  in  C ^ and  and  for  a non- trivial  solu- 

tion the  determinant  of  the  coefficients  must  vanish,  i.e.. 


(\m  + 2M.m)Tl(1)A;[1)  + X^A^1* 


SA™  + 


(Xm  + 2p-ln)Tl(2)A^2)  +Xm§A3(2) 


?a<2)  + V2)a3(2) 


Equation  (4.28)  is  a complex  algebraic  equation,  both  the  real  and  imaginary 

parts  of  which  must  vanish  simultaneously.  Solutions,  i.e.,  values  of  w 

and  § satisfying  (4.24)  and  (4.28)  are  found  numerically  in  a manner  similar 

to  that  discussed  in  the  paragraph  following  Eq.  (3.18). 

Calculations  have  been  performed  for  a set  of  material  parameters  cor- 

4 

responding  to  a glass  fiber  reinforced  phenolic  resin  , the  relevant  constants 
of  which  are 

pm  = 0.00013  lb-sec2/in4  , Ef  = 12.4  X 106  lb/in2, 

pf  = 0.00026  lb-sec2/in4  , |if  = 10.2  X 106  lb/in2, 

Xm  = 0.86  X 106  lb/in2,  p,1"  = 0.37  X 106  lb/in2,  (4.29) 

for  a fiber  diameter  of  .01  in.  for  the  volume  percentage  of  reinforcement  of 

5.67%,  which  corresponds  to  s = .04  inches.  The  results  of  the  calculations 
are  plotted  in  Fig. 4,  which  indicates  the  existence  of  an  upper  (optical 
type)  surface  wave  branch  in  addition  to  the  lower  (acoustic  type)  branch. 

As  already  noted,  we  do  not  believe  that  the  upper  surface  wave  branch 
actually  exists,  but  that  its  existence  is  a consequence  of  the  reduced 
coupling  in  the  simplified  model.  A similar  analysis  has  been  performed  for 
straight-crested  surface  waves  propagating  normal  to  the  direction  of  the 
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fiber  reinforcement  and  the  results  of  calculations  based  on  this  analysis 
are  plotted  in  Fig. 5.  Again  the  results  indicate  the  existence  of  an  upper 
as  well  as  a lower  branch  and  we  do  not  believe  that  the  upper  branch 
actually  exists  for  the  aforementioned  reasons.  In  Figs. 4 and  5 we  have 
drawn  vertical  lines  which  correspond  to  a wavelength  five  times  the  spacing 
of  the  fiber  reinforcement.  We  do  not  believe  the  curves  to  be  valid  much 
beyond  these  vertical  lines  because  of  the  nature  of  the  model  of  the  com- 
posite we  have  employed,  and  we  draw  them  considerably  beyond  their  range 
of  validity  simply  to  indicate  the  calculated  behavior.  The  important  curves 
in  Figs. 4 and  5 are  the  lower  acoustic  type  branches,  both  of  which  are 
drawn  to  a larger  scale  in  Fig. 6.  Note  the  difference  in  dispersion  for  the 
two  directions  of  propagation  considered.  This  very  precise  dispersion  property 
of  surface  waves  could  well  be  used  as  a means  of  nondestructively  evaluating 
the  distribution  of  the  fiber  reinforcement  in  and  the  integrity  of  the 
bonding  to  the  matrix. 
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APPENDIX 

In  this  Appendix  we  present  the  solution  to  the  second  problem  treated 
in  Sec. 2,  i.e.,  the  one  associated  with  Fig. 2,  using  the  reduced  equations 
presented  in  Sec. 4 and  omitting  the  influence  of  gravity.  In  this  case  the 
differential  equations  (4.15)  and  (4.16)  with  L = 3 and  the  constitutive 
equations  (4.6)  with  L = M = 3 and  (4.9)^,  all  in  the  absence  of  any  x^-, 
x^-dependence  and  inertial  terms  respectively  replace  Eqs.  (2.12),  (2.13), 

(2.2)  and  (2.4).  Equations  (2.27),  (2.28),  (2.30)  and  (2.32)  remain  unchanged 

and  (2.26)  and  (2.29)  are  replaced  by 


u3(1)  = u3(2)  , «“>  - U3  at  x3  = 0 , 


(1)  (2)  m 

T + T = T 

33  33  133 


at  x3  = 0 


The  solution  functions  given  in  (2.14)  and  (2.33)  are  replaced  by 

(1)  005  3 -a,X3 

U3  = C1  + C2X3  + C36  + C46  ’ U3  = C5  + C6X3  » 


“f  ' C1  + C2*3  + 


r,  2 (n*+^mn„  ox3  . -“3, 

L1  - * - -- ' j c3e  + c4e  > - 


where 


,3  m „ m (2).  . (2)  ,.m  _ m,,^ 

a = [a2( A + 2m.  + e )/E  (A  + 2m  )] 


(Al) 

(A2) 


(A3) 


(A  4 ) 


Substituting  from  (A3)  into  (Al),  (A2),  (2.30)  and  (2.32),  we  obtain 

P b 


- p 3 

c.  = — 2 . c = c*3p  , c = — 

, 9 > -1 


1 .m  m’  2 *0  '32  cosh  oil  ’ 

A + 2m 


A2 . 


The  actual  stresses  in  the  fibers  and  t™  in  the  matrix  at  the  junction  are 

JJ  JJ 


given  by 


. m „ m_ 


33 


33 


, m , _ m , _ 

X + 2[i  + E 


f (2)  \ + 2^1 

(2)  L cosh  oil  ' ’ 


(pV/A) (Am  + 2um) 
r + 2^+E(2) 


h _ 1 1 

L cosh  oil J ’ 


(A7) 


where 


pf  = p A/Af  , E(2)  = Ef Af/A  , 

o o ’ ’ 


(A8) 


and  p^  denotes  the  actual  stress  in  the  fibers  before  they  enter  the  matrix. 
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Schematic  Diagram  of  Loaded  Fiber  Reinforced  Composite 
Supported  on  a Rigid  Base 

Schematic  Diagram  of  Loaded  Fiber  Reinforced  Composite  with 
Reinforcement  Terminating  Uniformly  in  Matrix  Some  Distance 
Before  Support 

Schematic  Diagram  Showing  Surface  Wave  Propagating  Axong 
a Free  Surface  of  a Fiber  Reinforced  Composite  and  the 
Associated  Coordinate  System 

Dispersion  Curves  for  Surface  Waves  Propagating  in  the 
Direction  of  the  Fiber  Reinforcement 

Dispersion  Curves  for  Surface  Waves  Propagating  Normal 
to  the  Direction  of  the  Fiber  Reinforcement 

Acoustic  Type  Dispersion  Curves  for  Surface  Waves 
Propagating  Both  in  and  Normal  to  the  Direction  of  the 
Fiber  Reinforcement 
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